The European Physical Journal E

Healing of the gastrointestinal (GI) tract following surgical interventions is governed by coupled mechanobiological processes and may lead to severe complications such as perforation or stricture formation. A first comprehensive mechanobiological modeling framework of GI healing is herein proposed. The framework addresses two clinically relevant applications of bioprinting: endoscopic treatment of GI tissue defects and bioprinting-assisted surgical anastomosis. Tissue remodeling is based upon a hyperelastic material model with fiber dispersion coupling the spatiotemporal reaction diffusion dynamics of fibroblasts and cytokines with collagen deposition, reorientation, and plastic deformation. An active stress formulation drives fibroblasts and myofibroblast contraction while a multiplicative decomposition approach is adopted for collagen remodeling as a permanent change. The governing equations are discretised in time with an implicit Euler scheme, finite element in space and solved by a staggered algorithm implemented in the open-source platform FEniCS. Numerical simulations show that the initial collagen fibers content and orientation strongly influence fibroblast infiltration and scar contracture formation. The study quantifies the advantage of having a collagen bioprinted architecture aligned with tissue fibers, to reduce structural remodeling distortion. The model also predicts a strong correlation between collagen deposition and luminal burst pressure during the recovery of mechanical integrity in GI anastomosis.

High-speed atomic force microscopy (HS-AFM) experiments provide unique insights into nanoscale biological processes by directly visualizing the surface topography dynamics of biomolecules, such as proteins, over extended periods of time and generating large imaging data sets. However, the lack of computational methods for automated data analysis has thus far limited the full utilization of these data for the quantitative understanding of single-molecule dynamic. Here, we present the AFM Analysis Arena (A3) within BioAFMviewer as a software platform for automated, high-throughput quantitative analysis of AFM topographic imaging data within a user-friendly interactive graphical interface, implementing a convenient workflow from raw HS-AFM (AFM) imaging to statistical analysis. In several applications to HS-AFM imaging data, we demonstrate methods of automatized molecule detection and dynamic tracking to highlight the available capacities, aiming to transform HS-AFM measurements from observations towards an analytical tool in nanoscale biology.

Cardiomyocyte contraction is essential for the pumping action of the heart and deteriorates after myocardial damage, either as a consequence of irreversible cardiomyocyte injury or stiffening of the extracellular matrix, a process described as fibrosis. Cell geometry and substrate stiffness not only influence sarcomere architecture and contractility, they also determine how much work the cardiomyocytes can transfer to their environment. Human induced pluripotent stem cell (hiPSC)-derived cardiomyocytes provide a model system to examine these effects in a controlled manner, enabling us to mimic key features of the cardiac environment in vitro. Here we show that geometric confinement and substrate stiffness influence different aspects of cardiomyocyte force generation. By combining traction force microscopy and live-cell imaging of hiPSC-derived cardiomyocytes on soft and patterned substrates, we find that geometric confinement leads to longer resting sarcomeres and more pronounced sarcomere shortening on substrates with both physiological (10 kPa) and fibrotic (30 kPa) stiffness. Interestingly, the substrate stiffness itself has little effect on resting sarcomere length, but influences the peak contractile stress. On 10 kPa substrates, confined cells generate higher peak contractile stresses than unconfined cells, whereas this difference is not observed at 30 kPa. Confinement also results in faster mechanical relaxation at 10 kPa, with no detectable difference at 30 kPa. The combination of defined cell geometry and physiological stiffness leads to longer resting sarcomere lengths, more pronounced sarcomere shortening, higher peak contractile stress, and faster relaxation, features associated with a more mature cardiomyocyte phenotype. These findings show that the mechanical consequences of cell geometry depend on substrate stiffness and, therefore, both aspects have to be considered when employing hiPSC-derived cardiomyocytes as a model system in mechanobiology research, disease modeling and drug testing.

Cellulose has long been viewed as the major load-bearing component in plant primary cell walls. Recent studies, however, suggest that pectin plays a significant role in cell wall mechanics despite a kPa-scaled modulus. Here we quantify the mechanical properties of centimeter-sized internodal cell walls of Chara corallina following treatments to remove calcium crosslinks or digest pectin. We also characterize the flow of water and gas through the cell wall as a function of pressure. We found that removal of calcium crosslinks did not affect wall strength and stress-relaxation behavior, yet it decreased the elastic modulus and increased permeability to water. In contrast, pectin removal had a dramatic impact on all measurements. Samples from which pectin had been removed exhibited lower strength and stiffness and had greater stress relaxation under constant strain. Pectin removal increased the permeability of the cell wall to water above that observed when calcium was removed and reduced the pressure at which gas "tunnels" through the cell wall ~ 2-fold. Our findings support the idea that pectin locally restricts the movement of cellulose microfibrils, with implications for understanding functional properties of primary cell walls in vascular plants and developing high-performance double-network materials.

Biofilm mechanics depend strongly on growth history, hydration, nutrient availability, and interfacial attachment, yet conventional rheology requires harvesting and transfer that can perturb the as-grown structure. Here, we introduce bioRheoFlow, a 3D-printed rheometer-mounted platform for cultivation and mechanical characterization without sample transfer. The platform combines evaporation compensation with static or axial nutrient flow, and benchmark measurements validate its rheometric performance. Applied to Pseudomonas putida and Pantoea sp., bioRheoFlow reveals that biofilms cultivated in situ under axial flow exhibit substantially higher elastic moduli than conventionally grown and scraped samples. Despite similar linear-viscoelastic stress ranges, the species show distinct nonlinear responses: P. putida exhibits a lower crossover stress and more abrupt failure, whereas Pantoea sp. sustains larger stresses and deformations. Creep suggests contributions from bulk deformation and biofilm-substrate interactions. bioRheoFlow provides a controlled framework for probing growth-conditioned mechanics in living biofilms.

Rule-based languages have proven effective for modeling systems of interacting structured entities as typically encountered in chemistry and molecular biology. We present PyKappa, a rule-based modeling package written in Python whose interpreted nature enables interactive simulation and analysis, including by agentic AI. The package seeks to broaden the base of developers by utilizing a widely known programming language and serves as an easy-to-deploy teaching tool. Using PyKappa, we conduct a case study of phase separation.

Known for their rapid water surface locomotion, whirligig beetles (Gyrinidae) are also fast divers. They can transition from the surface to underwater swimming in under 100 ms. Diving requires them to generate sufficient hydrodynamic torque to break the water surface and reorient the body. However, the mechanism by which whirligigs create the necessary pitch-down torque is unknown. Here, we show that their diving is achieved by bending the abdomen. The inertial hydrodynamic pressure on the flexed abdomen generates the torque necessary for surface breaking and body rotation. The torque is maximized at an abdominal angle of 45{degrees}, whereas for effective diving, whirligigs display a preferential abdominal angle of around 28{degrees}, closely aligning with the angle that maximizes the torque-to-drag ratio. This abdominal angle is optimized by balancing torque generation with the drag penalty. Our findings inspire a mechanism for surface-to-underwater transition, offering design principles for multimodal robots capable of locomotion across the air-water interface.

Adherent cells must continuously adapt to rapid environmental fluctuations to preserve mechanical integrity. However, conventional theoretical frameworks predominantly rely on quasi-static assumptions, thereby limiting their ability to capture the transient dynamics of cellular responses to high-rate perturbations, such as acute osmotic shocks. To address this limitation, we developed a biophysical model incorporating a Hill-type law that couples active cytoskeletal mechanics to cell-edge velocity, explaining the counterintuitive experimental observation that traction forces transiently decrease during rapid hypotonic swelling despite an increase in cell size. We further construct a kinematic-geometric phase diagram that connects classical quasi-static theories with our dynamic model. Moreover, we show that traction-force dynamics are sensitive to the loading rate and amplitude of osmotic shock. In addition, stiffness-associated simulations and supporting perturbation evidence suggest a cellular stiffness-associated recovery trend in which stiffer cells exhibit faster volume recovery driven by stronger recoil of hydrostatic pressure.

Centriolar satellites are membrane-less organelles that store, modify and traffic proteins, in particular those of the centrosome and the primary cilium, and their deregulation is linked to ciliopathies, neurodevelopmental and neurodegenerative disorders. Satellite number, size, intensity and spatial distribution vary with cell type and cell state, yet satellites are usually quantified by a single readout, the pericentrosomal fluorescence intensity in maximum-projected images, which cannot distinguish changes in satellite number, size and position. Here we present SatelliteQ, an open-source Fiji and Python workflow for object-based quantification of satellites in 2D, 3D and time-lapse images. SatelliteQ segments satellites, detects centrosomes or primary cilia in the same cells, and reports satellite number, size, intensity, distance to the centrosome, shape and spatial distribution per satellite and per cell, with interactive graph and spatial viewers. Its measurements correlate closely with Imaris-based quantification; absolute values depend on segmentation settings and analysis mode, whereas the biological differences tested here were preserved across matched settings. Using SatelliteQ, we resolved distinct spatial signatures of satellite dispersal induced by microtubule depolymerization and by disruption of weak hydrophobic interactions, quantified differences between cell types and the scaling of satellite number with cell size, captured PCM1 granule dynamics in living cells, and compared satellite organization between ciliated and unciliated cells. Depletion of the ciliopathy-associated satellite protein CEP290 altered satellite abundance and pericentrosomal organization, while patient-derived cells carrying CEP290 mutations displayed divergent satellite phenotypes that were not explained by ciliation alone. SatelliteQ thus establishes satellite properties as a multidimensional cellular phenotype and provides an accessible framework for exploring satellites as an organelle-level readout in disease.

A cell's ability to maintain its volume is a prerequisite for its essential physiological functions. Changes in extracellular osmolarity induce the flux of water and solutes across the cell and nuclear membrane and, therefore, alters cell and nuclear volumes. Here, we investigated the equilibrium and dynamic volumetric responses of cells and nuclei in adherent human hepatoma (HepG2) cell monolayers using live three-dimensional imaging. Through image segmentation, 3D reconstruction, and volume rendering, precise cell and nuclear volumes were quantified. We found that during hypoosmotic stress, cells swell and lose osmolytes through active volume regulation, and nuclei undergo greater relative expansion than cells. During hyperosmotic stress, we found that both cells and nuclei shrink rapidly, develop irregular morphologies and undergo passive volume regulation. Furthermore, we explored the osmotic responses in these cells by restoring isosmotic conditions following osmotic stress and found that nuclei returned to volumes greater than their initial isosmotic volume, whereas cytoplasmic volume returned to a value below its initial isosmotic volume following swelling. By fitting coupled transport models water and solute permeability estimates for both cells and nuclei were made from their dynamic responses to different anisosmotic media with and without the presence of permeating solutes dimethyl sulfoxide and ethylene glycol. Together, these quantitative measurements of cell, cytoplasm, and nuclear volumetric responses and associated transport parameters have important implications for understanding cell physiology and the mechanism of osmoregulation.

Despite lacking a dedicated structural function, brain tissue is exposed to severe mechanical loading during accidents and trauma, with potentially permanent or lethal consequences. Therefore, characterizing its mechanical behavior assumes particular importance for designing protective equipment, planning surgeries, and training neurosurgeons, tasks that often rely on physical brain replicas made from phantom materials. Although the mechanics of brain tissue has been extensively studied, the selection of suitable phantom materials still relies on simplified metrics that fail to capture the highly nonlinear and rate-dependent behavior of the brain. In this work, we characterize the mechanical response at low-to-medium strain rates of common synthetic brain phantom materials, including Sylgard 527, Ecoflex Gel, and a composite hydrogel, and establish rigorous comparisons with porcine brain tissue as a benchmark. To this end, we employ a multi-step stress relaxation protocol under both uniaxial extension and compression across two different strain rates, introduce thermodynamic metrics, and perform comparisons based on rigorous statistical analysis. In particular, we isolate the thermodynamic equilibrium stress-strain curves, representing the baseline material's response, from rate-dependent dissipative inelastic effects, relatable to viscoelastic and poroelastic mechanisms. Our results provide a thorough comparison, revealing critical limitations in existing simulant materials and highlighting significant space for developing improved brain surrogates. The comprehensive dataset collected can also support the calibration of constitutive models, enhancing the predictive accuracy of numerical simulations involving brain tissue and its surrogates. More broadly, this study highlights the need to move beyond simplistic, scalar-metric-based comparisons of mechanical responses across materials.

Biofilm formation--the establishment of cellular communities within a self-secreted polymer matrix--is a behavior exhibited by almost all species of bacteria. Bacterial biofilms are well described as viscoelastic materials, with properties that grant physical advantages, such as elastic resilience and viscous adaptability. Additionally, many biofilm-forming species produce a subpopulation of highly resilient spores. In the model biofilm- and spore-forming species Bacillus subtilis, sporulation and matrix production are regulated by a common gene pathway, provoking the question of how sporulation influences biofilm viscoelastic properties. Investigating this point will facilitate the management of both salutary and deleterious biofilms, especially through manipulating their establishment and dispersion. Here, we investigate the interplay of sporulation and viscoelastic properties using rheological measurements of B. subtilis biofilms with varied matrix and spore production. We find that, to a significantly greater extent than matrix production, biofilm physical development is strongly altered by activation of sporulation through the spo0A regulation pathway. Spore-deficient strains initially establish physically robust biofilms with linear viscoelastic properties and recovery behavior distinct from spore-producing biofilms. However, spore-deficient biofilms fail to maintain their physical properties over time. In contrast, spore-producing wild-type and matrix-mutant biofilms exhibit more stable physical properties. Our results demonstrate that cell-level gene regulation is associated with macroscopic mechanical transitions in cellular communities, with sporulation granting not only a cellular-level survival advantage, but also greater community-level physical stability.

Tropomyosins play a central role in the regulation of the actin cytoskeleton by assembling along filaments into rope-like, head-to-tail polymers ("clusters") along filaments that regulate the binding and activity of numerous actin-binding proteins. Importantly, in cells, tropomyosins do not decorate actin filaments uniformly; some filaments are devoid of tropomyosins, whereas others recruit specific isoforms. To explore the diversity of tropomyosin-isoform binding kinetics and identify factors that control their assembly, we designed tropomyosin constructs, bearing a natural N-terminal acetylation and targeted cysteine mutations for fluorescence labelling. We analysed the binding of different tropomyosin isoforms to individual actin filaments reconstituted inside a microfluidic chamber. We show that, despite significant differences in spontaneous nucleation and elongation rates, all tested tropomyosin isoforms can coassemble within the same clusters. In addition, we identified tropomodulin and leiomodin, two pointed-end-binding proteins, to efficiently nucleate the formation of tropomyosin clusters on actin filaments, for all tested tropomyosin isoforms. Notably, while tropomodulin nucleates tropomyosin clusters from the pointed end, leiomodin does so along the entire filament.

We investigate the dynamics of the spread of the infection due to a phage- during the range expansion of an E coli colony. We find a rich panoply of spatio-temporal dynamics at large colony-level length scales, arising due to the interplay of mechanical forces, growth and viral infection dynamics at smaller cellular length scales. Strikingly, the phages "surf" the front of the growing bacterial colony by hitchhiking on E coli that are advected due to the growth of bacteria in the colony, resulting in an anisotropic spread of the phages. We identify microscopic processes - of the phage release during lysis (occurring on short millisecond timescales at length scales close to that of a single bacterium) and local nematic alignment of the rod-like E coli bacteria (occurring on timescales comparable to the bacterial growth rate and the length scales of a few bacterial cells) - that enhance the advective effects driving the hitchhiking behavior. Combining our experiments with a mathematical model, we explain the multiple spatio-temporal dynamical regimes - from the coexistence of uninfected bacteria, resistant (lysogenic) cells and phages, to fixed points where the entire population turns resistant or remains uninfected. Our work broaches a new frontier in the investigation of eco-evolutionary processes arising from the mechanical and ecological forces due to the interactions between viruses and their hosts.

Centrosome positioning is governed by mechanical interactions between the actin and microtubule cytoskeletons, yet the contribution of different actin architectures and their attachment to the cell membrane remains unclear. We use numerical simulations to study a centrosomal array of microtubules interacting with cortical, bulk, and centrosome-associated actin networks, including passive cross-linkers, motor proteins, and ERM-mediated coupling of actin to the cell membrane. We find that bulk and centrosome-associated actin alone are insufficient for robust centrosome centering, whereas sufficiently dense cortical actin promotes centering through steric interactions with growing microtubules. ERM-mediated attachment of cortical actin to the cell membrane enhances this effect, allowing stable centering at lower actin densities. In a two-centrosome system, microtubule interactions generate a spindle-like structure, while the organization of actin and its membrane attachment through ERM regulate centrosome separation. Antagonistic microtubule motors further limit centrosome separation and establish a finite spindle length. These results demonstrate that ERM-mediated coupling of actin to the cell membrane provides an important mechanical link between the actin cortex and microtubule-driven centrosome positioning, influencing both centrosome centering and spindle geometry.

Multi-cellular pattern formation is a key process in organism development, regeneration, and disease. Such multi-cellular patterns include, for instance, cell density and gene expression patterns. Here, we develop a cell-track-based method that quantifies how much individual cellular processes contribute to the emergence of large-scale patterns. These cellular processes include changes in cellular expression levels, cell motion, cell divisions, and apoptosis. We validate our method using synthetic data, before demonstrating it on two sets of experimental data. Our method is purely geometric and does not require any mechanistic assumptions about the underlying pattern-forming processes, making it very generic and robust. It thus enables unbiased, data-driven study of the mechanisms driving pattern formation across a broad class of systems.

Clathrin-mediated endocytosis (CME) is the primary pathway for internalization of extracellular and membrane cargo in eukaryotic cells. It is characterized by a patch of the plasma membrane invaginating and pinching off to generate a cytoplasmic vesicle. Recently, experimental evidence has begun to accumulate in support of proteins assembling into liquid-like condensates at nascent CME sites, but potential functional consequences are not understood. Given that condensate capillary forces can influence membrane bending, we constructed complementary analytical and numerical model systems to probe the mechanical influence of a protein droplet on endocytosis. Model evaluations and analysis reveal that a droplet can regulate the progression of endocytosis by stalling invagination until the droplet grows past a threshold size. After the droplet-mediated brake is released, the droplet provides a force that can drive the membrane bud toward scission. Thus, we propose a ``mechanical timer'' mechanism as a function for protein phase separation in cellular membrane trafficking.

Nanoplastic interactions with cell membranes can influence particle uptake, translocation, and biological effects, making their molecular description important for assessing the environmental and biological consequences of plastic pollution. However, molecular simulations generally represent nanoplastics using monodisperse polymer chains, whereas nanoplastics formed by environmental degradation are polydisperse. Here, we construct 4 nm polyethylene terephthalate (PET), polyethylene (PE), and polystyrene (PS) nanoparticles with lognormal molecular-weight distributions reported for degraded polymers, and with melt-like chain packing and entanglement. Using Martini coarse-grained molecular dynamics, we simulate their interactions for 10 microseconds with three lipid bilayers of increasing complexity. PE and PS insert rapidly and remain largely intact, whereas PET remains surface-associated and progressively releases individual chains. Fluid, polyunsaturated lipid environments undergo greater deformation and promote substantially more chain release, up to 56% of the chains for PET and 43% for PS. Release begins with the shortest chains and progressively extends to longer chains, while semicrystalline PE remains largely intact. In contrast, a size-matched monodisperse PS particle releases no chains and minimally perturbs the membrane. These results identify the molecular-weight distribution as a key determinant of nanoplastic-membrane interactions, and indicate that uniform model particles underestimate the release of polymer chains into biological membranes.

Oncogene amplification on circular extrachromosomal DNA (ecDNA) has been linked to poor prognosis and higher treatment resistance in multiple types of human cancer. ecDNA are mobile genetic elements lacking centromeres that are partitioned unevenly into daughter cells at mitosis. While random segregation of ecDNA contributes to gene-copy-number heterogeneity among tumour cells, how ecDNA contribute to phenotypic heterogeneity, via mRNA and proteins, and how this affects targeted therapy outcomes is still not understood. In fact, cancer cell populations with ecDNA show remarkable heterogeneity in mRNA concentrations, which is not explained by gene-copy-number heterogeneity alone. In this paper, we use a stochastic model of gene transcription to mathematically model mRNA production and decay in cancer cells with ecDNA, extending a method that has been used to study phenotypic heterogeneity in eurokaryotic and prokaryotic cells. Using our model, we explore the space of possible mRNA distributions from a given ecDNA distribution, under different model parameter regimes. Our model is able to produce realistic mRNA distributions and we use it to estimate transcriptional parameters from single-cell measurements of ecDNA and mRNA counts in cell-line data. We discover a remarkable similarity in estimated parameters across cell lines and amplification types, indicating similar processes underlying transcription in ecDNA and chromosomal amplifications. We also model genotoxic therapy and use our model to predict the mRNA composition of persisting cells, showing the usefulness of such a modelling framework to understanding the origins of treatment resistance in ecDNA-containing cancers.

Trigger waves are positive feedback-mediated traveling waves that propagate without slowing down or diminishing in amplitude, enabling biochemical coordination over distances where diffusion alone would be insufficient. Apoptosis can spread as trigger waves of caspase activity, as demonstrated in Xenopus eggs and egg extracts. We measured apoptotic wave speeds in extracts at two cell cycle states and found that waves propagated more slowly in M-phase arrested cytostatic factor (CSF) extracts than in interphase extracts (21 versus 27 m min-1). Luther's equation, v=2{surd}D{kappa}, predicts wave speed from the diffusion coefficient D and the positive-feedback doubling rate {kappa}. However, D and {kappa} were comparable between the two extracts and could not explain the speed difference. This discrepancy highlighted a key limitation of the classical equation: it does not account for an activation threshold in trigger wave initiation, and in fact the threshold required to trigger apoptosis was found to differ between cell cycle states. We therefore extended Luther's equation to include such a threshold, obtaining v{propto}{surd}(AD/Uc ), where A is the reaction rate and Uc is the activation threshold. Perturbing apoptotic sensitivity confirmed the predicted threshold-speed relationship. These findings establish the activation threshold as a key determinant of trigger wave speed and provide a broadly applicable scaling law for predicting the spread of switch-like biochemical signals.